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Abstract 



We consider perturbative solutions to the classical field equations coming 
from a quadratic gravitational lagrangian in four dimensions. We study the 
charged, spherically symmetric black hole and explicitly give corrections up 
to third order (in the coupling constant [j multiplying the R^uR^'^ term) to 
the Reissner-Nordstrom hole metric. We consider the thermodynamics of 
such black holes, in particular, we compute explicitly its temperature and 
entropy-area relation which deviates from the usual S = A/A expression. 
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I. INTRODUCTION 



Quadratic or fourth-order gravity appears to be a good candidate to represent the low 
energy hmit of the, yet unknown, quantum theory of gravity. In fact, it arises as the 
necessary counterterms one encounters when trying to renormahze semiclassical theory at 
one-loop level and as the low-energy limit of string theory (ignoring the dilaton and 
antisymmetric fields) [^f. Consequently, it is interesting to study the classical solutions of 
such effective theory. In Ref. we have developed a method to solve the field equations 
of the quadratic gravitational theory in four dimensions coupled to matter (see also Refs. 
^ where a closely related approach was independently derived from a formal perturbation 
theory). The quadratic terms are written as a function of the matter stress tensor and 
its derivatives in such a way to have, order by order, Einsteinian field equations with an 
effective T^u as source. By successive perturbations around a solution to Einstein Gravity, 
which for us represent the zeroth order, one can build up approximate solutions. In Ref. , 
we applied this perturbative procedure to find first order solutions in the coupling constant 
f3 for the charged black hole (where in this particular case a contributions vanish). 

We consider the following Lagrangian formulation of quadratic theories 

I = Ia + I^ = -^ f d^x^l ~2A + R 

lOTT J L 

+aR^ + PR^.R^"' + IQircJj , (1) 

where we have set G = 1 and c = 1 for simplicity. 

For the perturbative approach to properly work we consider only small curvatures, such 
that 

a\R\ < 1 , \pRf,^\ < 1 . (2) 

The coupling constants a and (3 are expected to be of the order of the Planck length 
squared. They already have upper bounds of the atomic size order from observations in the 
solar system, binary pulsars and cosmology. 

It is worth to stress that the coupling constants a and (3 must fulfill the, so called, 
no-tachyon constraints 

3a + /? > , /3 < , (3) 

that can be deduced upon linearization and asking for a real mass for both the scalar field 
related to R and the spin- two field '^p^lu related to R^i, (see Ref. [0 for further details). 

The field equations derived by extremizing the action I can be rewritten to the n-th 
order approximation as 

where the zeroth-order corresponds to the ordinary Einstein equations and 
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H^, = -2R,^, + 2g^,UR - -g^,R^ + 2RR^, , (5) 

and 

I^i, = -2R'^ .^^ + aR^^ + -g^^UR 

+2R°Rau — -gfj.uRaisR"'^ ■ (6) 

In Section II we explicitly compute the effective stress tensor, metric coefficients and 
horizon radius to third order in the coupling constant (3. In Section III we use this results 
to analyze the thermodynamic properties of charged black holes in the quadratic theory of 
gravity. We then brieffy discuss the validity of the laws of thermodynamics for black holes 
in theories of gravitation different from Einstein's general relativity. Finally, we discuss 
the relevance of the higher order /^-corrections in the computations of the extreme black 
hole case, and in the "stabilization" of the quadratic solution around the general relativistic 
metric. 



II. CHARGED BLACK HOLE METRIC 



We shall now study spherically symmetric solutions to the quadratic field equations 
representing charged black holes. Its metric can be written (in the Schwarzschild gauge) as 

ds"^ = gttir)df + grrir)dr^ + r'^d^f , (7) 

where dVt^ = d'd'^ + siv? 'dd({P' . The starting point is the General Relativistic solution, i.e. 
the Reissner-Nordstrom metric, —gtt{f) = grr{f)~^ — ~ 2M/r + Q^/r^). 
The non-vanishing components of the energy-momentum tensor are 

Tl = T:.=-Tt = -T^ = -^^. (8) 

The exact metric, solution of the Einstein equations with the effective source, Eq. (P, 
can be written as [g| [in the Schwarzschild gauge, Eq. (0)] 

g;,\r) = 1 + — / f^Tf 'dr , 9 

r r Joo 



gu{r) = -grr{r)-^ exp jsyr £ (T; - Tlr^fg„{f)df^ . (10) 

We use Eq. (|^) as the zeroth order effective T^^/ to compute the higher order (Eq. 
(D) we use a program of analytic manipulation within "Mathematica" |^ to find up to third 
order corrections to the energy-momentum tensor 
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45 - 220— + 280— + 135^ - 345^^ + 106^ 
-T^TTT^ 75600 - 635040— + 1774080^ + 430080-^ - 2432640^^ 
-1646400— + 847908-^ + 3385200 ^ - 2323965^^ + 525476-^ 

/yO ry\) ryK) J 

+0{P') , (11) 

* 87rr4 27rr6 V r r2 / 

_^f225 - 1140^1 + 1400^f + 745^i; - 1755^^ + 528^? 

OVT"?^ \ ry ryZi ry /yrj ry^ 



^ ^ ' 176400 - 1542240— + 4294080— + 1123920-^ - 5993680- ^ 



'^J^Qyj-^lO \ /y ry'2i ry'2. ry'3i 



-3841600— + 2005316-V + 7772240 — ^ - 5069085^^ + 1074388-^ 

ryO /y^ ly^ ryU ry\) . 

. (12) 



Since the effective energy-momentum tensor keeps its property of being traceless, the 
other two non-vanishing components can be deduced from the two above, i.e. 

Tf ^ = Tf^ = -^{Tf * + Tf , (13) 

Replacing this into Eqs. (il)-(|in) we obtain the metric components 

9rr\r) = 1 - — + ^ - ^ (lO - 15- + 6^ 

+ T7^ 7560 - 31920— + 33600— + 17880-^ - 36855^^ + 9856-^ 

J^^J^^O \ ry ry^ ry^ ryO ly^ 

-V A 14414400 - 110270160— + 272912640— + 71431360-^ 

5005r* V n^l n^-l 

-342838496^^ - 219739520— + 104276432-^ + 404156480 — -f- 

ryO /yO ry^ ly^ 

-241626385^^ + 47273072-^ + 0{(3^) , (14) 

/yO fy\) J 

Wr) = -fl- — + ^U^f5-5- + ^ 

8^Q^ ( M m2 g2 Mg2 g4 

^ ^ ' 2520 - 8400— + 6720— + 3600-^ - 5355^^ + 952^ 



^Q^^6 \ /y ry'2i ry'2i /yS ry^ 

/33g2/ M M2 13504 g2 
f^-^ - 720 + 4656 9600^ — 

ryO \ ry ry A ^\ ry ^ 
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592864 MQ2 3403504 0^ 112064 M^Q^ 
H ^ + 6272— H + 



55 1155 11 r 



2295703 MQ^ _ 12352965^ 64M^ 

429 r5 1365 91 ^ ' 'J ^ ' ^ ^ 

Since the complete expression for ga is already rather complicated we do not pursue the 
computation of higher order corrections. 

The validity of this metric will be assured if the condition Eq. @ holds. In our case, 
this takes the form —PQ'^/r^ <^ 1. Thus, Eqs. (|T5|)-(|14|) will be a good approximation to 
the (even extremely) charged black hole solutions in quadratic theories if 



th > ^J-P , (16) 

where th is the radial coordinate of the event horizon (in Schwarzschild's gauge). 

The radial coordinate of the horizon can be computed from making vanish gtt{fH) given 
by Eq. (p!5|). To first order in j3 we have (for non-extreme black holes) 



I3Q'^ (5r2 - 3Q2) ^2g4 (2<^2br\ - 6515Q2r^ + 5095QV2 _ ^gg^g^ 

rjj = r+ + 



5r3 {rl - g2) i050r| (^ - Q^f 

750750ri_i {rl - Q^^ 

36916630QV'^ - 435320Q^°r^ - 5714693Q^V2 + 1442637Q^^ ) + C»(/3^) , (17) 



600600r^^ - 43507050QV1_2 + 75284175QV"' - 79088475gV^ 



\5 



where r+ = M + ^M'^ - 

We observe that the quadratic black hole shrinks with respect to the corresponding 
general relativistic one when we consider up to terms, but the next order already tends 
to stabilize the solution near the general relativistic value. In fact, we can expect r^^ > for 
a wide range of /?, thus ensuring the existence of charged black holes in quadratic theories 
of gravitation. 

The extreme black hole will be now reached with a maximal charge different from the 
general relativistic one 

«L, = MV?^+A^ + 0(/3'), (18) 
thus Th given by Eq. (|l3), will remain always bounded. In fact, 

= M + ^ + ^(/^') • (19) 

[Note that since gu ~ (r — r'jf^^)'^ for an extreme black hole to obtain r'jf^^ to order jS"' we 
need to know g^ to order /J^"]. 

The horizon area, A^, is given by 
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167r/32Q4 



5254 -r^^^ 



+ , . : . ' , . , (-287g^ + 1069QV^ - 1340gV^ + 600r 

^ [434532gi^ - 999705Q^V2 - 5007445Qi°r^ + 24633770gV^ 



375375rf (g2 _ 

43546450gV^ + 38687775gV^° - 17253525gVi2 + 3003000ri* ) + 0{P^) . (20) 



III. THERMODYNAMIC PROPERTIES OF THE PERTURBATIVE SOLUTION 

In the theory of black holes in classical general relativity a mathematical analogy was 
discovered between certain laws of black hole mechanics and the ordinary laws of thermo- 
dynamics. With the discovery of quantum particle creation near black holes, this analogy 
appeared to be something more than a simple mathematical analogy. 

Very recently, some aspects of this significative relation between black hole mechanics 
and thermodynamics have been studied in a more general context than in Einstein's general 
relativity. In particular, here, we want to know for our approximate solution what is the 
explicit form of the laws of black hole mechanics. 

As it is known the zeroth law of black hole mechanics asserts that the surface gravity k 
(where n is defined by ^^Vf,^" = kC,"" with is the null generator of horizon) is constant all 
over the event horizon of a stationary black hole. The proof of this law makes direct use 



of the specific form of Einstein field equations and matter dominant conditions |jTO|, and it 
was not extended to others theories of gravity. Anyway, it can be easily shown that the 
constancy of the surface gravity trivially holds in the spherically symmetric case and can be 
computed using Eq. (^) below. 

In fact, we can safely compute the Bekenstein-Hawking temperature from the surface 
gravity, k, since this equation was found euclideanizing the black hole metric without making 
any assumption on the field equations that it fulfilled H], 



rp ^ 9tt 

J-H 



27r An ^J-gttgr 



(21) 



r=rH 



where g'^^ = dgu/dr. 

Thus, in our approximation, for non-extreme black holes we have, 

5257rr|_^ {Q^ — r^) ^ 

- 1848000g^^ - 46018351687giV 

^ 30030007rri.7(r2 -g2)n ^ ^ + 

+85649730010g^V^ + 215820442547g^V^|. - 752950763928g^°r^ 
+751221815775gV^° - 212346839950gV^2 _ 88152761275gV^'^ 

+46904894700gV^^ - 126126000r^^ + C(l/r+)') + C(/3^) , (22) 
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Note that, since /3 < 0, the effect of the quadratic gravitational theory corrections at first 
order in (3 will be that of decreasing the black hole radiation temperature with respect to the 
general relativistic value (with the same M and Q), leaving thus out open the possibility of 
switching off black hole evaporation and leaving behind a charged remnant with a mass of 
the order of the Planck mass. Higher order corrections, however, alternate sign. To further 
discuss this possibility one should have an exact solution and stablish the values of (3 that 
allow switching off of the black hole temperature. From Eq. (^) and Eq. (^) we can check 
that TniQmax) = up to order (3. 

The question of the first law of black hole mechanics in a general theory of gravity was 
largely discussed in Refs. |TT| , p!2| and clearly holds in our case, since the quadratic theory 
is derived from the diffeomorphism invariant lagrangian (0). Inverting and integrating the 
fundamental relation 



dM 



OA 



dA 



dQ 



dQ = TndS + ^ndQ , (23) 



A 



in Ref. 0, we found, to first order in /5, an approximate expression for the entropy of 
the charged black hole, in which clearly appears not simply one quarter of the area of the 
horizon, but in general a more complicated function of the area, A and charge, Q. 



Recently, in Ref. [|12[, Wald has developed a mathematical rigorous method to derive the 
exact formal expression of the entropy of a stationary black hole in any general covariant 
theory of gravity. In fact, the entropy was found to be a local geometrical quantity expressed 
as an integral evaluated on an arbitrary cross-section of a killing horizon over the associated 
Noether charge. 

In Ref. Wald's Noether charge techniques have been explicitly applied to a lagrangian 
theory such as (|l|), which is clearly generally covariant (invariant under general diffeomor- 
phism transformations). Their general expression for the entropy can be written as 

S=^J^ d^xVh[l + 2aR + (3g^^R^u] , (24) 

where S is any arbitrary cross-section of the horizon, d'^xy/h is the intrinsic volume element 
and where g'^ = g^" — h^^ = {^^x^ + X^^^) is the metric in the subspace normal to the 
horizon, with the stationary killing field and ^ vector field orthogonal S satisfying 
X''Xp = and = 1 on the horizon. 

In fact, making use of expression (P^), the entropy for our approximate solution takes 
the following form 



S = l f dM^^gM^i^ + Wg''Rtt\ = 7 fl + 167r/3Tf *! _ x / dM^^g^^ 

4: Jr=r[] 4 L At—vh Jr=rH 



r=rH 



= ^ _ Stt^/?^ - 10247r^/53^(A - 47rg2)2 + iQSSAir'' (3^ {A - AnQ^' x 

(25A^ - mTrQ^A^ + 35067r2Q4^2 _ iii20n^Q^A + 124487r^g^) + C(/5^) , (25) 

where we have taken the bifurcation surface as the cross-section of the horizon S and where, 
in the charged spherically symmetric case, R = and g^l = 2/ g^ is the only non vanishing 
component of the metric in the subspace normal to the horizon. In addition, Eq. ([T2|) was 
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used for computing the term g^^Ru = SnT^^ *, and the area of the event horizon is given by 
Eq. di. 

Note that the knowledge of the metric to allowed us to compute the entropy up 

to 0{f3^). That means that from the expression of the Reissner-Nordstrom metric we can 
compute S to 0{P), and thus recover the results of Refs. [0,|14| without having to solve any 
further field equation for g^iy. The vanishing of the term proportional to jS"^ is a casual fact 
due to the particular dependence of T^^ * to (that order vanishes when evaluated on 

the horizon). 

In Ref. [^] it was noted that, since the action Eq. (|ip "is linear in a and in jS, the 
modifications to the entropy from higher curvature terms in the original action would only 
be linear in a and in The precise meaning of this statement could not appear clear to 
the reader at first sight. When we write explicitly the entropy as a function of the horizon 
area and charge, we find corrections coming from quadratic curvature terms in the general 
formula Eq. that are not linear, but of {n + l)-th order in /3, if the metric contains n-th 
order terms in this coupling constant. 

To continue our previous discussion on the validity of the laws of black hole mechanics 
in the case of our approximate solution, we consider the the second law of thermodynamics, 
which states that the entropy must be an always increasing function. The theorem that 
establishes that the area of the black hole horizon must be an always increasing quantity, 
was proved only for general relativity |jl5 . 



Although any proof of the second law of black hole mechanics in a general theory of 
gravity has not been given yet, it seems to exist a relation between its validity and the 
positive energy condition of the theory. In the case of lagrangian theory examined in this 
paper (Eq. ([l|)), the positive energy condition does not holds in general, as we can see 
from the linearized field equations given in Ref. |]^, where the presence of the spin two 
massive tensorial field ip"''' allows for negative values of the energy conditions on the energy- 
momentum tensor. 

As in general relativity, here, the third law of thermodynamics should hold in the form 
due Nerst: Loosely speaking, "it is very hard to reach = 0". The version due to Planck 
does not hold, since from Eqs. ([T8|)-(p!9|) we have 



S(Th = 0) = ttM^ -—(3 + 0{(3^IM^). (26) 

5 

This ground state is clearly degenerate due to the general relativistic term. The contribution 
of quadratic theories is (to first order) only to contribute with an universal constant (function 
of/3). 

With regards to the fourth law or the scaling laws fulfilled by the black hole critical 



exponents [16,1/], we expect them to hold also for our approximate solution. In partic- 



ular, the black hole critical exponents should be the same as in general relativity, since 
the Reissner-Nordstrom solution gives the leading behavior and since the validity of the 
universality hypothesis [jl7|. 



IV. DISCUSSION 

The aim of this work was two-fold: First to see if, in practice, our perturbative approach 
allowed to be carried out to orders higher than the first. This has been shown to be possible 
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explicitly up to third order. The second objective of this work was to see if from the higher 
orders one could guess some recursive formulae that would lead us to an exact solution. So 
far, we have not fulfilled this expectation, but may be the reader can do it. We can, however, 
try to say something about the /3-gravity (or /3-sector of the quadratic theories): Although 
the first order corrections seem to indicate that /3-gravity is weaker than General Relativity, 
we can see (for example in Qmax, fH and Th) that to higher order we have alternateness in 
the sign of the corrections (we needed to go up to 3rd. order in to check that). Using 
the Noether charge method we have been able to express the black hole entropy in terms 
of the horizon area and charge up to 4th. order in /?, and (in spite of the casual vanishing 
of the 0{l3^) term) we have shown that the dependence of 5* on /5 is not simply linear. We 
have also studied the extremely charged black hole. There, due to the particularity of this 
solution, to obtain r'jf^^ (and the other quantities evaluated on the extreme horizon such as 
Tfj and S) to the first order in /3 we had to compute g^^, to order P^. 

From the fact that our perturbative method can not go to the high curvature regime, 
unfortunately, we cannot say anything about the most interesting issue of the singularity at 
r = 0. 

We are presently dealing with the rotation effects, i.e. the f3 first-order corrections to the 
Kerr-Newman solution and the results will be published elsewhere |TH . 
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